The mechanism and high-free-energy transition state of lac repressor-lac operator interaction.

Significant, otherwise-unavailable information about mechanisms and transition states (TS) of protein folding and binding is obtained from solute effects on rate constants. Here we characterize TS for lac repressor(R)-lac operator(O) binding by analyzing effects of RO-stabilizing and RO-destabilizing solutes on association (ka) and dissociation (kd) rate constants. RO-destabilizing solutes (urea, KCl) reduce ka comparably (urea) or more than (KCl) they increase kd, demonstrating that they destabilize TS relative to reactants and RO, and that TS exhibits most of the Coulombic interactions between R and O. Strikingly, three solutes which stabilize RO by favoring burial/dehydration of amide oxygens and anionic phosphate oxygens all reduce kd without affecting ka significantly. The lack of stabilization of TS by these solutes indicates that O phosphates remain hydrated in TS and that TS preferentially buries aromatic carbons and amide nitrogens while leaving amide oxygens exposed. In our proposed mechanism, DNA-binding-domains (DBD) of R insert in major grooves of O pre-TS, forming most Coulombic interactions of RO and burying aromatic carbons. Nucleation of hinge helices creates TS, burying sidechain amide nitrogens. Post-TS, hinge helices assemble and the DBD-hinge helix-O-DNA module docks on core repressor, partially dehydrating phosphate oxygens and tightening all interfaces to form RO.

Fluorescence Resonance Energy Transfer Characterization of DNA Wrapping in Closed and Open Escherichia coli RNA Polymerase-λP(R) Promoter Complexes.

Initial recognition of promoter DNA by RNA polymerase (RNAP) is proposed to trigger a series of conformational changes beginning with bending and wrapping of the 40-50 bp of DNA immediately upstream of the -35 region. Kinetic studies demonstrated that the presence of upstream DNA facilitates bending and entry of the downstream duplex (to +20) into the active site cleft to form an advanced closed complex (CC), prior to melting of ∼13 bp (-11 to +2), including the transcription start site (+1). Atomic force microscopy and footprinting revealed that the stable open complex (OC) is also highly wrapped (-60 to +20). To test the proposed bent-wrapped model of duplex DNA in an advanced RNAP-λP(R) CC and compare wrapping in the CC and OC, we use fluorescence resonance energy transfer (FRET) between cyanine dyes at far-upstream (-100) and downstream (+14) positions of promoter DNA. Similarly large intrinsic FRET efficiencies are observed for the CC (0.30 ± 0.07) and the OC (0.32 ± 0.11) for both probe orientations. Fluorescence enhancements at +14 are observed in the single-dye-labeled CC and OC. These results demonstrate that upstream DNA is extensively wrapped and the start site region is bent into the cleft in the advanced CC, reducing the distance between positions -100 and +14 on promoter DNA from >300 to <100 Å. The proximity of upstream DNA to the downstream cleft in the advanced CC is consistent with the proposed mechanism for facilitation of OC formation by upstream DNA.

E. coli RNA Polymerase Determinants of Open Complex Lifetime and Structure.

In transcription initiation by Escherichia coli RNA polymerase (RNAP), initial binding to promoter DNA triggers large conformational changes, bending downstream duplex DNA into the RNAP cleft and opening 13bp to form a short-lived open intermediate (I2). Subsequent conformational changes increase lifetimes of λPR and T7A1 open complexes (OCs) by >10(5)-fold and >10(2)-fold, respectively. OC lifetime is a target for regulation. To characterize late conformational changes, we determine effects on OC dissociation kinetics of deletions in RNAP mobile elements σ(70) region 1.1 (σ1.1), ß' jaw and ß' sequence insertion 3 (SI3). In very stable OC formed by the wild type WT RNAP with λPR (RPO) and by Δσ1.1 RNAP with λPR or T7A1, we conclude that downstream duplex DNA is bound to the jaw in an assembly with SI3, and bases -4 to +2 of the nontemplate strand discriminator region are stably bound in a positively charged track in the cleft. We deduce that polyanionic σ1.1 destabilizes OC by competing for binding sites in the cleft and on the jaw with the polyanionic discriminator strand and downstream duplex, respectively. Examples of σ1.1-destabilized OC are the final T7A1 OC and the λPR I3 intermediate OC. Deleting σ1.1 and either ß' jaw or SI3 equalizes OC lifetimes for λPR and T7A1. DNA closing rates are similar for both promoters and all RNAP variants. We conclude that late conformational changes that stabilize OC, like early ones that bend the duplex into the cleft, are primary targets of regulation, while the intrinsic DNA opening/closing step is not.

Quantifying additive interactions of the osmolyte proline with individual functional groups of proteins: comparisons with urea and glycine betaine, interpretation of m-values.

To quantify interactions of the osmolyte l-proline with protein functional groups and predict their effects on protein processes, we use vapor pressure osmometry to determine chemical potential derivatives dµ2/dm3 = µ23, quantifying the preferential interactions of proline (component 3) with 21 solutes (component 2) selected to display different combinations of aliphatic or aromatic C, amide, carboxylate, phosphate or hydroxyl O, and amide or cationic N surface. Solubility data yield µ23 values for four less-soluble solutes. Values of µ23 are dissected using an ASA-based analysis to test the hypothesis of additivity and obtain α-values (proline interaction potentials) for these eight surface types and three inorganic ions. Values of µ23 predicted from these α-values agree with the experiment, demonstrating additivity. Molecular interpretation of α-values using the solute partitioning model yields partition coefficients (Kp) quantifying the local accumulation or exclusion of proline in the hydration water of each functional group. Interactions of proline with native protein surfaces and effects of proline on protein unfolding are predicted from α-values and ASA information and compared with experimental data, with results for glycine betaine and urea, and with predictions from transfer free energy analysis. We conclude that proline stabilizes proteins because of its unfavorable interactions with (exclusion from) amide oxygens and aliphatic hydrocarbon surfaces exposed in unfolding and that proline is an effective in vivo osmolyte because of the osmolality increase resulting from its unfavorable interactions with anionic (carboxylate and phosphate) and amide oxygens and aliphatic hydrocarbon groups on the surface of cytoplasmic proteins and nucleic acids.

Quantifying why urea is a protein denaturant, whereas glycine betaine is a protein stabilizer.

To explain the large, opposite effects of urea and glycine betaine (GB) on stability of folded proteins and protein complexes, we quantify and interpret preferential interactions of urea with 45 model compounds displaying protein functional groups and compare with a previous analysis of GB. This information is needed to use urea as a probe of coupled folding in protein processes and to tune molecular dynamics force fields. Preferential interactions between urea and model compounds relative to their interactions with water are determined by osmometry or solubility and dissected using a unique coarse-grained analysis to obtain interaction potentials quantifying the interaction of urea with each significant type of protein surface (aliphatic, aromatic hydrocarbon (C); polar and charged N and O). Microscopic local-bulk partition coefficients K(p) for the accumulation or exclusion of urea in the water of hydration of these surfaces relative to bulk water are obtained. K(p) values reveal that urea accumulates moderately at amide O and weakly at aliphatic C, whereas GB is excluded from both. These results provide both thermodynamic and molecular explanations for the opposite effects of urea and glycine betaine on protein stability, as well as deductions about strengths of amide NH--amide O and amide NH--amide N hydrogen bonds relative to hydrogen bonds to water. Interestingly, urea, like GB, is moderately accumulated at aromatic C surface. Urea m-values for protein folding and other protein processes are quantitatively interpreted and predicted using these urea interaction potentials or K(p) values.

One-step DNA melting in the RNA polymerase cleft opens the initiation bubble to form an unstable open complex.

Though opening of the start site (+1) region of promoter DNA is required for transcription by RNA polymerase (RNAP), surprisingly little is known about how and when this occurs in the mechanism. Early events at the lambdaP(R) promoter load this region of duplex DNA into the active site cleft of Escherichia coli RNAP, forming the closed, permanganate-unreactive intermediate I(1). Conversion to the subsequent intermediate I(2) overcomes a large enthalpic barrier. Is I(2) open? Here we create a burst of I(2) by rapidly destabilizing open complexes (RP(o)) with 1.1 M NaCl. Fast footprinting reveals that thymines at positions from -11 to +2 in I(2) are permanganate-reactive, demonstrating that RNAP opens the entire initiation bubble in the cleft in a single step. Rates of decay of all observed thymine reactivities are the same as the I(2) to I(1) conversion rate determined by filter binding. In I(2), permanganate reactivity of the +1 thymine on the template (t) strand is the same as the RP(o) control, whereas nontemplate (nt) thymines are significantly less reactive than in RP(o). We propose that: (i) the +1(t) thymine is in the active site in I(2); (ii) conversion of I(2) to RP(o) repositions the nt strand in the cleft; and (iii) movements of the nt strand are coupled to the assembly and DNA binding of the downstream clamp and jaw that occurs after DNA opening and stabilizes RP(o). We hypothesize that unstable open intermediates at the lambdaP(R) promoter resemble the unstable, transcriptionally competent open complexes formed at ribosomal promoters.

Probing DNA binding, DNA opening, and assembly of a downstream clamp/jaw in Escherichia coli RNA polymerase-lambdaP(R) promoter complexes using salt and the physiological anion glutamate.

Transcription by all RNA polymerases (RNAPs) requires a series of large-scale conformational changes to form the transcriptionally competent open complex RP(o). At the lambdaP(R) promoter, Escherichia coli sigma(70) RNAP first forms a wrapped, closed 100 bp complex I(1). The subsequent step opens the entire DNA bubble, creating the relatively unstable (open) complex I(2). Additional conformational changes convert I(2) to the stable RP(o). Here we probe these events by dissecting the effects of Na(+) salts of Glu(-), F(-), and Cl(-) on each step in this critical process. Rapid mixing and nitrocellulose filter binding reveal that the binding constant for I(1) at 25 degrees C is approximately 30-fold larger in Glu(-) than in Cl(-) at the same Na(+) concentration, with the same log-log salt concentration dependence for both anions. In contrast, both the rate constant and equilibrium constant for DNA opening (I(1) to I(2)) are only weakly dependent on salt concentration, and the opening rate constant is insensitive to replacement of Cl(-) with Glu(-). These very small effects of salt concentration on a process (DNA opening) that is strongly dependent on salt concentration in solution may indicate that the backbones of both DNA strands interact with polymerase throughout the process and/or that compensation is present between ion uptake and release. Replacement of Cl(-) with Glu(-) or F(-) at 25 degrees C greatly increases the lifetime of RP(o) and greatly reduces its salt concentration dependence. By analogy to Hofmeister salt effects on protein folding, we propose that the excluded anions Glu(-) and F(-) drive the folding and assembly of the RNAP clamp/jaw domains in the conversion of I(2) to RP(o), while Cl(-) does not. Because the Hofmeister effect of Glu(-) or F(-) largely compensates for the destabilizing Coulombic effect of any salt on the binding of this assembly to downstream promoter DNA, RP(o) remains long-lived even at 0.5 M Na(+) in Glu(-) or F(-) salts. The observation that Esigma(70) RP(o) complexes are exceedingly long-lived at moderate to high Glu(-) concentrations argues that Esigma(70) RNAP does not dissociate from strong promoters in vivo when the cytoplasmic glutamate concentration increases during osmotic stress.

Interactions of the osmolyte glycine betaine with molecular surfaces in water: thermodynamics, structural interpretation, and prediction of m-values.

Noncovalent self-assembly of biopolymers is driven by molecular interactions between functional groups on complementary biopolymer surfaces, replacing interactions with water. Since individually these interactions are comparable in strength to interactions with water, they have been difficult to quantify. Solutes (osmolytes, denaturants) exert often large effects on these self-assembly interactions, determined in sign and magnitude by how well the solute competes with water to interact with the relevant biopolymer surfaces. Here, an osmometric method and a water-accessible surface area (ASA) analysis are developed to quantify and interpret the interactions of the remarkable osmolyte glycine betaine (GB) with molecular surfaces in water. We find that GB, lacking hydrogen bond donors, is unable to compete with water to interact with anionic and amide oxygens; this explains its effectiveness as an osmolyte in the Escherichia coli cytoplasm. GB competes effectively with water to interact with amide and cationic nitrogens (hydrogen bonding) and especially with aromatic hydrocarbon (cation-pi). The large stabilizing effect of GB on lac repressor-lac operator binding is predicted quantitatively from ASA information and shown to result largely from dehydration of anionic DNA phosphate oxygens in the protein-DNA interface. The incorporation of these results into theoretical and computational analyses will likely improve the ability to accurately model intra- and interprotein interactions. Additionally, these results pave the way for development of solutes as kinetic/mechanistic and thermodynamic probes of conformational changes and formation/disruption of molecular interfaces that occur in the steps of biomolecular self-assembly processes.

Late steps in the formation of E. coli RNA polymerase-lambda P R promoter open complexes: characterization of conformational changes by rapid [perturbant] upshift experiments.

The formation of the transcriptionally competent open complex (RP(o)) by Escherichia coli RNA polymerase at the lambda P(R) promoter involves at least three steps and two kinetically significant intermediates (I(1) and I(2)). Understanding the sequence of conformational changes (rearrangements in the jaws of RNA polymerase, DNA opening) that occur in the conversion of I(1) to RP(o) requires: (1) dissecting the rate constant k(d) for the dissociation of RP(o) into contributions from individual steps and (2) isolating and characterizing I(2). To deconvolute k(d), we develop experiments involving rapid upshifts to elevated concentrations of RP(o)-destabilizing solutes ("perturbants": urea and KCl) to create a burst in the population of I(2). At high concentrations of either perturbant, k(d) approaches the same [perturbant]-independent value, interpreted as the elementary rate constant k(-2) for I(2)-->I(1). The large effects of [urea] and [salt] on K(3) (the equilibrium constant for I(2) is in equilibrium with RP(o)) indicate that a large-scale folding transition in polymerase occurs and a new interface with the DNA forms late in the mechanism. We deduce that I(2) at the lambda P(R) promoter is always unstable relative to RP(o), even at 0 degrees C, explaining previous difficulties in detecting it by using temperature downshifts. The division of the large positive enthalpy change between the late steps of the mechanism suggests that an additional unstable intermediate (I(3)) may exist between I(2) and RP(o).

Solute probes of conformational changes in open complex (RPo) formation by Escherichia coli RNA polymerase at the lambdaPR promoter: evidence for unmasking of the active site in the isomerization step and for large-scale coupled folding in the subsequent conversion to RPo.

Transcription initiation is a multistep process involving a series of requisite conformational changes in RNA polymerase (R) and promoter DNA (P) that create the open complex (RP(o)). Here, we use the small solutes urea and glycine betaine (GB) to probe the extent and type of surface area changes in the formation of RP(o) between Esigma(70) RNA polymerase and lambdaP(R) promoter DNA. Effects of urea quantitatively reflect changes in amide surface and are particularly well-suited to detect coupled protein folding events. GB provides a qualitative probe for the exposure or burial of anionic surface. Kinetics of formation and dissociation of RP(o) reveal strikingly large effects of the solutes on the final steps of RP(o) formation: urea dramatically increases the dissociation rate constant k(d), whereas GB decreases the rate of dissociation. Formation of the first kinetically significant intermediate I(1) is disfavored in urea, and moderately favored by GB. GB slows the rate-determining step that converts I(1) to the second kinetically significant intermediate I(2); urea has no effect on this step. The most direct interpretation of these data is that recognition of promoter DNA in I(1) involves only limited conformational changes. Notably, the data support the following hypotheses: (1) the negatively charged N-terminal domain of sigma(70) remains bound in the "jaws" of polymerase in I(1); (2) the subsequent rate-determining isomerization step involves ejecting this domain from the jaws, thereby unmasking the active site; and (3) final conversion to RP(o) involves coupled folding of the mobile downstream clamp of polymerase.

The effects of upstream DNA on open complex formation by Escherichia coli RNA polymerase.

Binding of activators to upstream DNA sequences regulates transcription initiation by affecting the stability of the initial RNA polymerase (RNAP)-promoter complex and/or the rate of subsequent conformational changes required to form the open complex (RP(O)). Here we observe that the presence of nonspecific upstream DNA profoundly affects an early step in formation of the transcription bubble. Kinetic studies with the lambdaP(R) promoter and Escherichia coli RNAP reveal that the presence of DNA upstream of base pair -47 greatly increases the rate of forming RP(O), without significantly affecting its rate of dissociation. We find that this increase is largely due to an acceleration of the rate-limiting step (isomerization) in RP(O) formation, a step that occurs after polymerase binds. Footprinting experiments reveal striking structural differences downstream of the transcription start site (+1) in the first kinetically significant intermediate when upstream DNA is present. On the template strand, the DNase I downstream boundary of this early intermediate is +20 when upstream DNA is present but is shortened by approximately two helical turns when upstream DNA beyond -47 is removed. KMnO(4) footprinting reveals an identical initiation bubble (-11 to +2), but unusual reactivity of template strand upstream cytosines (-12, -14, and -15) on the truncated promoter. Based on this work, we propose that early wrapping interactions between upstream DNA and the polymerase exterior strongly affect the events that control entry and subsequent unwinding of the DNA start site in the jaws of polymerase.

The exclusion of glycine betaine from anionic biopolymer surface: why glycine betaine is an effective osmoprotectant but also a compatible solute.

Paradoxically, glycine betaine (N,N,N-trimethyl glycine; GB) in vivo is both an effective osmoprotectant (efficient at increasing cytoplasmic osmolality and growth rate) and a compatible solute (without deleterious effects on biopolymer function, including stability and activity). For GB to be an effective osmoprotectant but not greatly affect biopolymer stability, we predict that it must interact very differently with folded protein surface than with that exposed in unfolding. To test this hypothesis, we quantify the preferential interaction of GB with the relatively uncharged surface exposed in unfolding the marginally stable lacI helix-turn-helix (HTH) DNA binding domain using circular dichroism and with the more highly charged surfaces of folded hen egg white lysozyme (HEWL) and bovine serum albumin (BSA) using all-gravimetric vapor pressure osmometry (VPO) and compare these results with results of VPO studies (Hong et al. (2004), Biochemistry, 43, 14744-14758) of the interaction of GB with polyanionic duplex DNA. For these four biopolymer surfaces, we observe that the extent of exclusion of GB per unit of biopolymer surface area increases strongly with increasing fraction of anionic oxygen (protein carboxylate or DNA phosphate) surface. In addition, GB is somewhat more excluded from the surface exposed in unfolding the lacI HTH and from the folded surface of HEWL than expected from their small fraction of anionic surface, consistent with moderate exclusion of GB from polar amide surface, as predicted by the osmophobic model of protein stability (Bolen and Baskakov (2001) J. Mol. Biol. 310, 955-963). Strong exclusion of GB from anionic surface explains how it can be both an effective osmoprotectant and a compatible solute; analysis of this exclusion yields a lower bound on the hydration of anionic protein carboxylate surface of two layers of water (>or=0.22 H(2)O A(-)(2)).

Preferential interactions of glycine betaine and of urea with DNA: implications for DNA hydration and for effects of these solutes on DNA stability.

Interactions of the solutes glycine betaine (GB) and urea with mononucleosomal calf thymus DNA in aqueous salt solutions are characterized by vapor pressure osmometry (VPO). Analysis of osmolality as a function of solute and DNA concentration yields the effect of the solute on the chemical potential, mu(2), of the DNA. Although both GB and urea generally are nucleic acid denaturants and therefore must interact favorably with the nucleic acid surface exposed upon melting, VPO demonstrates that neither interacts favorably with duplex DNA. Addition of GB greatly increases mu(2) of DNA, indicating that the average local concentration of GB in the vicinity of the double helix is much less than its bulk concentration. By contrast, addition of urea has almost no effect on mu(2) of duplex DNA, indicating that the average local concentration of urea in the vicinity of duplex DNA is almost the same as in bulk solution. Qualitatively, we conclude that the nonuniform distribution of GB occurs primarily because duplex DNA and GB prefer to interact with water rather than with each other. Comparison with thermodynamic data for the interaction of GB with various protein surfaces (Felitsky et al., Biochemistry, 43, 14732-14743) shows that GB is excluded primarily from anionic DNA surface and that the hydration of anionic DNA phosphate oxygen surface (>or approximately 17 H(2)O per nucleotide or >or approximately 0.22 H(2)O A(-)(2)) involves at least two layers of water. From analysis of literature data for effects of urea and of GB on DNA melting, we propose that urea is an effective nonspecific nucleic acid denaturant because of its favorable interactions with the polar amide-like surface of G, C, and especially T or U bases exposed in denaturation, whereas GB is a specific GC denaturant because of its favorable interaction with G and/or C surface in the single-stranded state.

Preferential interactions in aqueous solutions of urea and KCl.

A quantitative characterization of the thermodynamic effects due to interactions of salt ions and urea in aqueous solution is needed for rigorous analyses of the effects of changing urea concentration on biopolymer processes in solutions that also contain salt. Therefore, we investigate preferential interactions in aqueous solutions containing KCl and urea by using vapor pressure osmometry (VPO) to measure osmolality as a function of the molality of urea (component 3) over the range 0.09

Kinetic studies and structural models of the association of E. coli sigma(70) RNA polymerase with the lambdaP(R) promoter: large scale conformational changes in forming the kinetically significant intermediates.

The kinetics of interaction of Esigma(70) RNA polymerase (R) with the lambdaP(R) promoter (P) were investigated by filter binding over a broad range of temperatures (7.3-42 degrees C) and concentrations of RNA polymerase (1-123 nM) in large excess over promoter DNA. Under all conditions examined, the kinetics of formation of competitor-resistant complexes (I(2), RP(o)) are single-exponential with first order rate constant beta(CR). Interpretation of the polymerase concentration dependence of beta(CR) in terms of the three step mechanism of open complex formation yields the equilibrium constant K(1) for formation of the first kinetically significant intermediate (I(1)) and the forward rate constant (k(2)) for the conformational change converting I(1) to the second kinetically significant intermediate I(2): R + P-->(K(1))<--I(1)(k(2))-->I(2). Use of rapid quench mixing allows K(1) and k(2) to be individually determined over the entire temperature range investigated, previously not possible at this promoter using manual mixing. Given the large (>60 bp) interface formed in I(1), its relatively small binding constant K(1) at 37 degrees C at this [salt] (approximately 6 x 10(6) M(-1)) strongly argues that binding free energy is used to drive large-scale structural changes in polymerase and/or promoter DNA or other coupled processes. Evidence for coupling of protein folding is provided by the large and negative activation heat capacity of k(a)[DeltaC(o,++)(a)= -1.5(+/-0.2)kcal K(-1)], now shown to originate directly from formation of I(1) [DeltaC(o)(1)= -1.4(+/-0.3)kcal K(-1)] rather than from the formation of I(2) as previously proposed. The isomerization I(1)-->I(2) exhibits relatively slow kinetics and has a very large temperature-independent Arrhenius activation energy [E(act)(2)= 34(+/-2)kcal]. This kinetic signature suggests that formation of the transition state (I(1)-I(2)++ involves large conformational changes dominated by changes in the exposure of polar and/or charged surface to water. Structural and biochemical data lead to the following hypotheses to interpret these results. We propose that formation of I(1) involves coupled folding of unstructured regions of polymerase (beta, beta' and sigma(70)) and bending of promoter DNA (in the -10 region). We propose that interactions with region 2 of sigma(70) and possibly domain 1 of beta induce a kink at the -11/-12 base pairs of the lambdaP(R) promoter which places the downstream DNA (-5 to +20) in the jaws of the beta and beta' subunits of polymerase in I(1). These early interactions of beta and beta' with the DNA downstream of position -5 trigger jaw closing (with coupled folding) and subsequent steps of DNA opening.